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Starting from the assumption that the present accelerated expansion of the Universe is driven 
by a chameleon scalar field, the function describing a direct coupling between the matter content 
of the Universe and the scalar field is derived from cosmological considerations. In progressing 
from cosmology towards astrophysical scales, this function is used to estimate the influence that 
such nonminimal coupling may have on the properties of usual (polytropic) stars whose matter 
can interact directly with the cosmological chameleon. It is shown that in order to obtain sizes 
and masses of compact configurations supported by the chameleon scalar field comparable to those 
of usual stars, it is necessary to assume that the pressure of the matter content of the Universe 
is substantially different from zero. The question of a possible manifestation of the effects of the 
f> | nonminimal coupling in laboratory experiments is also discussed. 
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PACS numbers: 95.36. + x, 95.35. + d, 04.40. - b 



I. INTRODUCTION 



The numerous attempts to give a theoretical description of the current accelerated expansion of the Universe are 
usually reduced to introducing an exotic form of matter called dark energy, which is believed to be responsible for such 
an acceleration. Whereas the true nature of dark energy is still unknown, different possibilities are being suggested to 
describe it. Perhaps one of the most promising approaches to addressing the origin of dark energy is a consideration 



of theories that include various fundamental fields Q1[H . Possibilities are also being considered .including some that 



involve modified (non-Einstein) gravity theories @,^] or models with extra space dimensions jlj, |(| . In any 



case, a 



description of the present acceleration of the Universe within the framework of the aforementioned approaches implies 
i ^ i , that one needs to choose the parameters entering these theories in such a way as to satisfy the current astronomical 
observations and laboratory experiments. 

One of such approaches is the so-called chameleon paradigm, the essence of which, briefly, is as follows 0, [1[ : Based 
on the assumption that scalar fields with a very small mass of order Hq ~ 10~ 33 eV do exist in the Universe, and they 
drive the current accelerated expansion, it is assumed that such fields may interact directly with ordinary matter and 
dark matter with gravitational strength. This leads to the fact that the scalar field acquires a mass which depends on 
^sjj j the local background matter density. For example, on Earth, where the density is high, the mass of the field becomes 
• . 10 30 times greater than the mass of this field on the cosmological background. 

An important ingredient of the chameleon mechanism is the presence of two functions of the scalar field ip: the 
potential energy V(tp) and the function f((p) describing the direct coupling between the matter and the scalar field. 
■ These two functions are arbitrary, and they are chosen in such a way as to satisfy various laboratory and cosmological 

^ ; tests 0-53- 

It is obvious that this arbitrariness in the choice of the mentioned field functions is due to lack of knowledge of the 
true nature of scalar fields filling the Universe (of course, assuming they do really exist in the present Universe and 
are responsible for its accelerated expansion). Perhaps the only thing one can do at this point is to try to find the 
restrictions imposed upon the form of these functions and values of the parameters appearing in them. Some work 
has already been done along these lines, particularly in the papers |UU - [l3j . where the authors choose a generalized 
expression for the nonminimal coupling between the scalar field and the matter having the form f{ip)L m , with L m 
being the Lagrangian of ordinary matter. In considering the evolution of the Universe they have used two different 
approaches: in Ref. the functions V(ip) and /(<£>) are taken to be arbitrary, and their parameters are chosen in 
such a way as to satisfy the current observational data. Another approach has been applied in [lj], HH , where the 
authors have initially selected some particular form of cosmological evolution and found the functions V(ip) and f(tp) 
corresponding to such evolution. 

On the other hand, the question may be asked: if the mechanism of direct coupling between a chameleon scalar 
field and ordinary matter is universal, i.e., it can be applied to any form of matter (be it ordinary matter or dark 
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matter), how will this affect the properties of various small-scale objects such as, for instance, usual stars? In Einstein 
gravity, compact configurations consisting of ordinary matter and scalar fields have been repeatedly considered earlier 
in the literature. The scalar fields may be self-interacting or not, and one can involve a fluid which interacts with the 
scalar field either only gravitationally [H], EH or through direct coupling as well [HI, E3] . In particular, in the recent 
papers [l8l - l20j the changes in the structure of compact objects (polytropic stars) under the influence of a chameleon 
field have been investigated. It was shown that, depending on the properties of arbitrarily chosen functions V(<p) 
and f(tp), the properties of these configurations, such as total mass, distribution of matter, and size, depend strongly 
on the surrounding scalar field. 

The purpose of the present paper is to demonstrate the possibility of determining the nonminimal coupling function 
f(ip) which satisfies the current astronomical observations. Subsequently, this function will be used in modeling 
compact objects in order to clarify the influence of the presence of the chameleon scalar field on their physical 
characteristics. As an example, a polytropic star embedded in the chameleon scalar field will be considered. 

The paper is organized as follows: In Sec. [II] the nonminimal coupling function f(ip) is obtained, in general form, 
from cosmological considerations. In the particular case of the ACDM model this function is derived in Sec. Ill Cl and. 
in Sec. IIII( is applied to the estimation of the influence of the nonminimal coupling upon the properties of usual stars 
whose matter can interact directly with a chameleon scalar field. In the Appendix, we consider the effects related to 
the nonminimal coupling in the sense of their possible manifestation in laboratory experiments. Finally, in Sec. IIVI 
we summarize and discuss the results obtained. 



II. / FROM COSMOLOGY 

As pointed out in the Introduction, here the choice of the nonminimal coupling function f(ip) will be done based on 
a consideration of the cosmological evolution of the Universe within the chameleon paradigm suggested in Ref. [TTI ] . 
To do this, we choose the Lagrangian of the system in the form 

L = --^-R+^d^d^-V^) + f(ip)L m . (1) 
l07r(_x z 

Here ip is the real scalar field with the potential V((f); A = ±1 corresponds to the usual or phantom scalar field, 
respectively; L m is the Lagrangian of ordinary matter (a perfect fluid). The case / = 1 corresponds to no direct 
coupling between the fluid and scalar field, but even in this case the two sources are still coupled via gravity. 

A. About the Lagrangian L m 

The Lagrangian L m can be chosen in a variety of ways. For example, in Ref. j2l| the choice L m — ~2e is used, and 
in Refs. [22|, |23|, L m — p (e and p are the energy density and the pressure of a fluid, respectively). By varying both 
these matter Lagrangians with respect to a metric, one obtains the same energy-momentum tensor of the perfect fluid 
in the conventional form. However, it can be shown that, for instance, for static configurations, these Lagrangians, 
being substituted in the general Lagrangian ([!}, will give different equations for an equilibrium configuration (the 
Tolman-Oppenheimer-Volkoff (TOV) equations). In the case of L m = p, the TOV equation will have the same form 
as when / = 1. On the other hand, the use of L m = —2e leads to the appearance of an extra term on the right-hand 
side of the TOV equation associated with the nonminimal coupling. 

Of course, in the limit / — > 1, for both choices of L m , the gravitating system of a scalar field plus ordinary matter 
reduces to the same uncoupled system in which the ordinary matter and the scalar field interact only gravitationally. 
But in the general case of / = f{(p), these two systems will not be equivalent, giving equilibrium configurations 
with different properties. Apparently, at the present time, it is difficult to make a motivated choice amongst these 
Lagrangians L mi or any other Lagrangians used in the literature (for the other possible Lagrangians, see, e.g., Ref. 24]). 
For this reason, a description of various systems with nonminimal coupling of the type f{(p)L m is made, in effect, 
using an ad hoc choice of L m . 

B. Deriving / in general form 

In this paper we choose L m = p which has been used in Refs. (l8l - [20j in modeling astrophysical objects. Here this 
choice of L rn will be used in studying the cosmological evolution of the present Universe. 
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Consider the flat model with the metric 



ds 2 = c 2 dt 2 -d\t)dV 



(2) 



where a(t) is the scale factor. Using the general Lagrangian (J]), the corresponding energy- momentum tensor is (details 
are given in the appendix of [Hj|) 



T? = / [(e + P )u lU k - 5$p] + A dap^ip - 5% 



d^ip - V(<p) 



(3) 



where e and p are the energy density and the pressure of the fluid, and u l is the 4- velocity. Using the metric ([2} and 
the energy-momentum tensor ([3]), the (q) component of the Einstein equations is then given by 
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A , 



(4) 



Here the dot denotes differentiation with respect to the cosmological time t. The equation for the scalar field <p follows 
from the Lagrangian ([T]) 
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and gives in the metric ((2J 
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One more equation can be obtained from the law of conservation of energy and momentum, T k k 
i = component of this equation gives 

d(fs) „a „. . df 
dt a dt 



(5) 



(6) 



0. Taking the 



(7) 



Keeping in mind that e and p are related by some equation of state, there are three unknown functions: a, ip, s, and 
three equations (H]), ©, and (J7J for them. 

Now we need to choose an equation of state for the matter content of the Universe. At the present time, it is 
believed that the main part of matter in the Universe (except dark energy) is in the form of dark matter. The most 
popular hypothesis in modeling dark matter is the assumption that its equation of state corresponds to cold (dust) 
matter Such matter is either pressureless or has a very small pressure (see, e.g., Ref. [Hj]). In hydrodynamical 
language, in cosmology, the relation between the pressure and energy density is usually taken in the simplest form 
as p — we, where w = const. Assuming w ~ at the present time, a description of the evolution of the Universe is 
carried out. 

Here we initially start with the assumption that the pressure of the matter filling the Universe is small compared 
with its energy density. To apply the results obtained from cosmology for modeling compact astrophysical objects, 
let us choose a polytropic equation of state which is widely used in astrophysics. This equation of state may be taken 
in the following parametric form: 



p = Kp 1+1 / n , e 



(8) 



Here the constant n is called the polytropic index and is related to the specific heat ratio 7 via n — 1/(7— 1), p 
is the rest-mass density, and K is an arbitrary constant whose value depends on the properties of the fluid under 
consideration. Substituting this equation of state in ([7]) , its general solution will be 



(1 + ap 1 ^ 



/a 3 ' 



(9) 



where po is an integration constant. The smallness of the pressure implies that the term in the denominator ap 1 ^ 71 <C 1. 
Then, neglecting this term, we obtain from Eq. Q the following approximate expression for the mass density: 



(10) 
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In the limit / — > 1, i.e., in the absence of the nonminimal coupling, this expression reduces to the usual dependence 
for the density of dust matter as a function of the scale factor. 

For the subsequent calculations, it is convenient to introduce new dimensionless variables 
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1/2 



¥>(*). 



(11) 



where B has the dimensions of inverse time. Using these variables and taking into account the equation of state 
and the expression (fT0|) , the system of Eqs. (UJ) and ([6]) takes the form 



»(-) + f + * 



(12) 
(13) 



Here the prime denotes differentiation with respect to the dimensionless time £; 7 = 1 + 1/n, V = V/(poc 2 ), a — 

Po/(poc 2 ) — Kp^ n /c 2 , where po,Po refer to the current values of the mass density and the pressure, respectively. 
Note that in our low-pressure approximation the parameter u«l. Also, in this approximation, the first term on the 
right-hand side of Eq. ([12")) always plays the role of dustlike matter, independently of the form of the function /. In 
turn, the latter, together with V , will describe an evolution of all other forms of matter, including dark energy. 

Next, using Eqs. (fT2j) and (p~3|) . let us find the form of the function / expressed in terms of the scale factor a. To 
do this, rewrite Eq. (TH?|) as 



Expressing from (|T2 



ld(cf>' 2 ) a' 
2 d£ + \ ( 



= 2A 




(14) 



a I a J 



and substituting this in fjl4[) . one can find the following dependence of / on a: 

f l -< = 6^ J a'a 3 ^ L 2 a" + 2aa' 2 - i - a 3 V^j d£. 



(15) 



(16) 



Then, if a and V are known as functions of the time £, it is possible to find the corresponding dependence / = /(£). 



C. Expression for / in the case of the ACDM model 

We now will demonstrate how to obtain an expression for / in the particular case of V = (a massless scalar field) , 
and when the ansatz for the scale factor a is chosen in the form of the well-known ACDM model [l| 




0(0 = 00 7^ Sinh ^f\/7T^ • (I 7 ) 



Here f2 m = 8itGpo/3Hq, Ha = c 2 A/3H§, where Ho — H(t — to) and ao = a(t — to) are the present values of 
the Hubble parameter and the scale factor, respectively. This expression for a smoothly interpolates between a 
matter-dominated Universe in the past (a oc t 2 / 3 ) and accelerated expansion in the future (a cx exp \/A/3 cb). 
Setting ao = 1 at the present epoch and substituting (fTTJ) in (|16p . one can find 



r 1-7 




(18) 
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In turn, we have from (|15l) 

(Here we have set A = +1 in order to ensure that the functions / and a obtained below will be real.) Note here 
that by substituting the above derivative of the scalar field into Eq. (|12p and multiplying it by Q m , we get the known 
equation for the ACDM model. That is, in the case under consideration, the gradient of the scalar field plays the role 
of the A-term. 

Substituting the expression for £ from ([!§)) in Eqs. (fTTj) and (fT5|) . we find the final formulas for / and a as functions 
of the scalar field <j> 




Using (|20|) . one can find the current value of the scalar field, <fio = <p(t = to). To do this, setting a = cio = 1 and 
choosing, for definiteness, £7a = 0.7 and tt m = 0.3, from the solution of the transcendental equation (f20|) . we have 

4> « 1.975, £ ~ 0.915. (22) 

This £q corresponds to to w 12.95 Gyrs. 



III. CHAMELEON STAR 



We now will proceed to a consideration of a compact configuration (a star) consisting of ordinary (polytropic) 
matter with the equation of state ([8]) and of a scalar field nonminimally interacting with this matter. In Ref. 18] this 
type of gravitating configuration of a polytropic fluid nonminimally coupled to a scalar field was called a "chameleon 
star." 

As a nonminimal coupling function, we will use the expression (|21j) obtained in the preceding section. We proceed 
from the assumption that the cosmological chameleon scalar field penetrates into the star and interacts with the star's 
matter through the same universal coupling function / from (|21[) . That is, we assume that the matter of the star, 
which is embedded in the external, cosmological, homogeneous scalar field, feels its presence not only gravitationally, 
but also through the nonminimal coupling. Then, the structure of the star is also defined by the distribution of the 
scalar field, and correspondingly by /, inside the star. 



A. General equations for hydrostatic equilibrium configurations 

We choose the static spherically symmetric metric in the form 

ds 2 =e"c 2 dt 2 -e x dr 2 -r 2 dn 2 , (23) 

where v and A are functions of the radial coordinate r, and dfl 2 is the metric on the unit two-sphere. The (§) and (}) 
components of the Einstein equations for the metric (|23|) and the energy-momentum tensor (J3j> are then given by 

OS = -e- (± ~ + i = ^ (/. + \ e" V + V) , (24) 





+ 1 


8ttG 


7) 




c 4 




+ 1 


8ttG 


7) 




c 4 



where the prime denotes differentiation with respect to 7'. (As in the previous section, we have again set A = +1.) 
Here the energy density s and the pressure p refer now to the matter of the star (but not to the cosmological matter), 
but the scalar field, in the spirit of the chameleon paradigm, remains the same as in cosmology. Our purpose will be 
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to clarify the influence that this cosmological scalar field may have on the inner structure and characteristics of such 
compact objects as usual stars. 

Equations (|24|) and ([25]) must be supplemented by an equation for the scalar field, which follows from ([5]) and has 
the following form in the metric (f2"3"|): 



2 1 
r 2 



dV 

dip 



dip 



One more equation follows from the law of conservation of energy and momentum, T^ k 
component of this equation gives 

dp 1 ^ dv 
dr 2 dr 

To perform the subsequent calculations, it is convenient to rewrite the equation of state ((5J in terms of the new 
variable 9 related to the density p at the given point and the central density p c by (see, e.g., Ref. [26|]) 



(26) 

0. Taking the i = 1 
(27) 



p = p c 6 n ° => p = Kp~<° = Kp l+1,n ° = Kp L c 
Notice that 9 must be "normalized" to unity when r = 0. Using (|28p in Eq. (|27p leads to 

d9 „ % df 



l+l/n a «n s +l 



2a s {n s + 1) 



<//• 



(28) 



(29) 



with a s = KpiJ n ° J (? — Pc/{PcC 2 ) and p c is the pressure of the fluid at the center of the configuration. The index 
s ("star") on a, n, and 7 is introduced here to distinguish these parameters from those used in Sec. [II] in describing 
the matter content of the Universe. In both cases, the relativistic parameters a and a s characterize the ratio of the 
pressure at a given point (in the case of a star, at the center of the configuration), or at a given instant of time (in 
the case of cosmology, at the present time) to the energy density. Since the expression for the coupling function / 
from (|2ip has been obtained in the low-pressure approximation, i.e., is valid only when a <C 1, therefore, it is natural 
that in the case of a star this expression is also applicable only when a s <C 1. 
Equation (|2U1) may be integrated to give e v in terms of 9 



I 



f 



2(n a +l) 



(30) 



where e Vc is the value of e" at the center of the configuration where 9 C — 1. The integration constant v c corresponds 
to the value of v at the center of the configuration. It is determined by requiring e" = 1 at infinity, i.e., that the 
space-time be asymptotically flat. 
We introduce a new function M(r) 



1 



2GM(r) 



M(r) is the mass of the configuration in the range < r < ri, where T\ is the boundary of the fluid where 
Also, we introduce dimensionless variables 

1/2 



x = Ar, v{x) 



A 3 M(r) 
4irp c 



<P(x) 
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<p(r), where A 



AiiGpc 



a s (n s + l)c 2 



_<7 s (n s + l)c 4 _ 

where A has the dimensions of inverse length. With this, we can rewrite Eqs. (|24[) and ([25]) in the form 

2 ~> 



dv 
dx 

, I — 2a s (n s + l)v/x d9 



f9 n s 



1 r 
2 



1 - 2a s {n s + 1) 



V 



1 



dr 



fe n - (1 - a. 



(31) 
= 0. 

(32) 

(33) 
(34) 



where V = V/(p c c 2 ) is the dimensionless potential energy of the field. Next, using ([2^]). one can rewrite Eq. (|2"r?)) as 
follows: 

d 2 cj) | (2 a s (n s + 1) \d0 

dx 



dx 2 



I 



1 + aJ 



1 f dv 



1 — 2a s (n s + l)v/x x \dx x 



dx 



1 - 2a s {n s + 1)- 



dV 



-idf 



(35) 
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Thus, the static configuration under consideration is described by the three equations (|33|) - (|35j) . 



B. The particular choice of / in the form (|21[) 

We now apply Eqs. (I33|) -(l35 j) for describing the inner structure of our configuration. In Ref. [llj], this system was 
used in studying similar configurations (chameleon stars) when the coupling function / was chosen arbitrarily. It 
was shown that the presence of the nonminimal interaction leads to substantial changes both in the radial matter 
distribution of the stars and in the star's total mass. 

To demonstrate the influence of the nonminimal coupling derived in Sec. |TT]on the properties of a compact config- 
uration, here we will consider a case when the potential energy is absent, i.e., V(</>) = 0. As mentioned above, we 
assume that the scalar field, being initially homogeneous (i.e., cosmological), in the presence of star's matter, begins 
to interact with the matter both gravitationally and through the nonminimal coupling. Our purpose will be to study 
the possible effects of this interaction. 

Note that, in the absence of the nonminimal coupling, the choice V(cj>) = corresponds to the case of a massless 
scalar field. However, one can see from the right-hand side of Eq. (I35|) that the presence of the nonminimal coupling 
leads to the presence of an additional term of the form a s 9 ns+1 df / 'dtp. 

Choosing / in the form (p?Tj) , Eqs. (|3"3")) - (|3"51) are to be solved for given a s , n s , and also cosmological parameters a, 
n, £Ia, and VL m subject to the boundary conditions in the vicinity of the center of the configuration x = 0, 

9^9 c + ^x 2 , v^v 3 x 3 , fa^fa + ^-x 2 , (36) 

where <f> c corresponds to the central value of the scalar field <fi, 9 C is normalized to be unity at the center, 9 C = 9(x = 
0) = 1. The values of the coefficients fa, i>3, and #2 are defined from Eqs. (|55|) . ([55]). and (|34l) . respectively, as 

fa = -\ a >{j£) ' "3 = /cA #2«-y(l + 4a s ), 

where the index c corresponds to central values. 

Since the cosmological scalar field, fat), ancl that used in modeling a star, fax), are measured in different units, 
then, to find the relation between them, we compare their definitions given by expressions (fTT|) and (|32|) . respectively. 
Thus we have 

fat) = ^2a s (n s + l)fax). (37) 
Substituting this expression in ([21]) , we finally obtain 

x -2(n+l) 



o 

- - Ill 



a(n + 1) 



sinli \/^cr s (n s + 1) fax) I . (38) 



As a result, the static configuration under consideration is described by the system of Eqs. (|33|) - (l35t together 
with the boundary conditions ([31))) and the function / from (f3"5]) . These equations, however, cannot be integrated 
analytically, and have to be solved numerically. The solutions were started near the vicinity of the center of the 
configuration (i.e., near x « 0) and solved out to a point x = x%, where the function 9 became zero. Since our 
spherical fluid configuration was taken to be embedded in an external, homogeneously distributed scalar field, we also 
needed to require that at the boundary, x = x%, the value of the varying scalar field from the inside of the star, fax), 
matched the exterior value fa x t- In this case, of course, the gradient of the scalar field, generally speaking, although 
small, is not strictly equal to zero. It requires continuing the solutions outside the fluid. Formally, the gradient goes 
to zero only asymptotically, i.e., as x — > 00. However, the numerical calculations indicate that for small values of the 
parameters a and a s used in this paper, the gradient of the scalar field is virtually equal to zero along the radius 
of the configuration. Then the form of the external scalar field is determined by matching it to the interior solution 
at the surface of the star. The junction conditions, which follow from the requirement that the field equations be 
satisfied on the boundary of the star, require that fax) be continuous. In our case of small scalar-field gradients, this 
can be done with sufficient accuracy by equating the interior and exterior values of the scalar field at x = Xi, i.e., by 
requiring that fax±) ss C xt- 

In turn, for a complete description of the system under consideration, it is also necessary to solve external Einstein's 
equations for the metric functions v and A from (|23p . This can be easily done in the manner described in Ref. |2fj| . 
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However, for the problems being considered in this paper, knowledge of the behavior of the exterior solutions is not 
necessarily needed. That is why, in order not to encumber the paper, we do not show this procedure here. 

Thus, the procedure for finding solutions is as follows: for given values of er, <r s , f2 m , n, and n s , we choose the 
central value of the scalar field, <j> c , such that goes to zero and 4>{x\) ~ cx t at some finite value of a; = x\. The 
fact that the configuration is embedded in the external cosmological scalar field 4>q implies that ox t must be equal 
to 4>q/ \j2cr s (n s + 1), where <fio is taken from (|22|) [the coefficient y/2<r s (n s + 1) appears because of the relation ([37]) ]. 
This allows one to determine definitely the value of <p c for such configurations. 

One can see from the expression (|38[) that the function / can be positive, negative, or imaginary depending on the 
value of n. Physically reasonable values of / are either positive or negative. However, the numerical calculations show 
that the parameters used in the paper give finite size solutions only for positive values of this function. To provide 
the positivity of /, it is necessary to take n either even or equal to 



21 



2m + 1 



(39) 



where m and I are integers. Proceeding from this choice, we get the results described below. 

In order not to encumber the paper, in further calculations the polytropic indexes n and n s are taken to be equal to 
each other. The calculations indicate that even if this is not the case, the qualitative behavior of the solutions remains 
the same. Keeping this in mind, we present here the results of numerical calculations with the following values of 
the parameters n and n s : n = n s — 6/5; 2; 8/3 (for convenience, we hereafter drop the index s on n s ). These values 
cover the range 1.2 < n < 2.67, which is often used in modeling usual stars. The parameters fl\ and fi TO are chosen 
to be JIa = 0.7, Cl m = 0.3, as before. In the approximation used here, the relativistic parameters u, <r s are taken to 
be much less than unity: a < 0.1, a s < 0.05. 

With this choice of the parameters, the numerical calculations indicate that the scalar field 4>{x) and the coupling 
function / are virtually constant along the radius of the configuration under consideration. This allows us to rewrite 
the system with / from (|38p, using the approximation (4'»i0 and f(x) ~ / c = const, in the following form: 



dv 
dx 

1 - 2a s (n + l)v/x dd 



1 



dx 



f c x 2 n , 

X 3 \fc0 n (1 - <Ts 



1 dv 

x 2 dx 



(40) 
(41) 



Then, introducing the new dimensionless radial coordinate x = xj \fY c and mass function v = u/v/c> we are able to get 
rid of f c in these equations. Thus, mathematically this system of equations becomes equivalent to that of describing 
usual polytropic stars with an equation of state in the form © and in the absence of a scalar field. Such configurations 
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were investigated in detail in the work of Tooper [27| who found regular, relativistic starlike configurations. These 
starlike solutions of [27| were found for the values of the parameters < a s < 0.75 and 1 < n < 3. 

In our case we are dealing with approximately the same range of n, but with small values of a s . Correspondingly, 
the size and the mass of our configurations are l/^/Jc times the values of those considered in [27j in the absence of a 
scalar field. That is, in our approximation the value of f c is the primary determining parameter of the problem. 

It is obvious that f c must depend on the values of three parameters: n, a, and a s . However, the numerical 
calculations show that, for fixed n and varying a and <r s , the argument of the hyperbolic sine in (|38[) is a slowly 
varying function of a s . That is why, to a good approximation, this argument may be regarded as a constant, that 
allows to rewrite the function / c (er) in the following approximate form: 



/cO) « /„ 



1) 



with /„ 




sinh 



-2(r, 



a s (n+ 1) 4>(x) 



const, 



(42) 



where /„ is approximately constant along the radius of the configuration for any fixed n. For the values of n used 
here, the coefficient /„ has the following numerical values: f n =6/5 ~ 0.342, f n —2 ~ 0.174, and f n =8/3 ~ 0.101. The 
numerical calculations indicate that, using these values, the approximation for f c (cr) in the form (|42[) is applicable 
with good accuracy when a < 0.1 and a s < 0.05. 

Figure Q] shows the function / c (ct) from (f42|) . One can see from this figure that the value of / c (c), being always much 
less than unity, grows with increasing a for all n. Since the size and the mass of the configurations under consideration 
are inversely proportional to the square root of / c (c), then its small values for any a and n correspond to large sizes 
and masses compared with configurations without a scalar field considered in Ref. |27j . The numerical results obtained 
in (27j are written in terms of dimensionless variables in the form (|11[) for arbitrary values of the equation-of-state 
parameters K, n, and the central density p c . Correspondingly, the concrete values of physical characteristics (masses 
and sizes) can be obtained in dimensional form by giving the corresponding values of the mentioned parameters. 
But in any case, the size and the mass of configurations with a scalar field considered in the present paper increase 
catastrophically in the limit of very small a, corresponding to dustlike cosmological matter. 



IV. SUMMARY AND DISCUSSION 



It is now believed that various scalar fields may exist in the present Universe. Based on the assumption that they 
can interact directly (nonminimally) with ordinary matter and dark matter with gravitational strength, this paper 
studies the possible influence that such fields may have on the properties of compact gravitating configurations (stars) 
consisting of ordinary matter. The choice of the nonminimal coupling function / is made from a consideration of the 
evolution of the Universe within the framework of chameleon cosmology. Then the total Lagrangian initially contains 
a term describing a nonminimal coupling between the scalar field tp and the matter in the general form f(ip)L m , with 
L m being the Lagrangian of ordinary matter [see Eq. ((TJ]. By choosing the arbitrary function f(ip), one might try to 
describe the evolution of the present Universe. 

In the present paper the choice of this function is based on two propositions: (i) The pressure of the matter filling the 
Universe is small compared with its energy density, (ii) The evolution of the cosmological scale factor a is determined 
by the observational data. 

To model the matter content of the Universe (namely, ordinary matter plus dark matter), we choose a polytropic 
equation of state in the form flSJ. This choice is caused by our wish to use the results obtained from cosmology in 
modeling compact objects (see below) . For this equation of state, in the low-pressure approximation, it is possible to 
find a general expression for / as a function of a in the form ([T5|. Next, using the well-known ACDM ansatz for the 
scale factor, we find / as a function of the scalar field in the final form (|2ip . 

Our next step is to consider the influence that the obtained nonminimal coupling may have on the properties of 
usual stars. While so doing, we start from the assumption that a star's matter is described by the same equation 
of state (jSJ) and also can interact directly with a cosmological scalar field. One can consider such a configuration as 
embedded in the scalar field which penetrates into the star and interacts with its matter through the function /. 

Depending on a concrete form of /, such a direct coupling results in more or less significant changes in the structure 
and properties of a star (see Refs. [l8l - [20l p. In our low-pressure approximation, when the pressures both of the 
cosmological matter and of the star's fluid are much less than their energy densities, the function /, as well as the 
scalar field, may be regarded as practically constant along the radius of the star: f(x) « f c = const. This allows 
one to rewrite the system of Eqs. (|33p - (135l) , describing a static configuration, in the form which is mathematically 
identical to that without a scalar field [see Eqs. (|40p and (|4ip]. This corresponds to the fact that the size and the mass 
of our configurations are \j\pf c times the values obtained in Ref. [27j . where configurations with the same equation 
of state, but without a scalar field, have been considered. But since in our case the values of f c are always much 
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less than unity, the size and the mass of configurations with the chameleon scalar field are always much greater than 
those of configurations without a scalar field. It is seen from Fig. Q] that this effect is especially strong for small values 
of the cosmological parameter a which correspond to the vanishing pressure of the matter filling the Universe (dark 
matter plus ordinary matter) . This gives rise to the fact that the size and the mass of compact configurations with 
the chameleon scalar field increase catastrophically compared with configurations without a scalar field. 

Does this unusual behavior of the characteristics of the compact configurations considered here inevitably imply 
that such chameleon cosmology is nonviable? The answer to this question might be given in the following way: As one 
can see from Fig. [TJ f c increases with increasing <j, and correspondingly the size and the mass of our configurations 
decrease. One might expect that the subsequent growth of <r, when the low-pressure approximation is not already 
valid, and the function / is not constant throughout the star, will provide the possibility of obtaining physically 
reasonable sizes and masses of configurations supported by the chameleon scalar field. At least the previous studies 
with an arbitrary choice of / performed in Refs. |18l - t20j | indicate that one can obtain sizes and masses comparable to 
those of usual stars without a scalar field. However, such a problem with / derived from cosmological considerations 
must be studied separately beyond the framework of the low-pressure approximation used in the present paper. 

Thus, within the framework of chameleon cosmology considered here, polytropic stars supported by the cosmological 
chameleon scalar field will have physically reasonable sizes and masses only if we suppose that the pressure of the 
matter filling the Universe is nonvanishing. However, in the case of relatively small pressures considered here, this 
does not lead to any significant changes from the point of view of conventional cosmology. Indeed, it is seen from the 
cosmological Einstein equation (|12[) that its right-hand side is identical in form to the source of matter in the form 
of dust plus the terms containing the scalar field which describe dark energy, regardless of its form. In the particular 
case of the ACDM ansatz for the scale factor considered here, the massless scalar field gave the contribution which is 
equivalent to the A-term. 

In the case when pressures both of the matter content of the Universe and of star's matter are not small compared 
with their energy densities, one would expect that the situation will change. On the one hand, one can expect to get 
more reasonable sizes and masses of compact configurations. On the other hand, from the point of view of cosmology, 
matter sources of the Einstein equations would be expected generally to have a different form from that of usually 
used in modeling the evolution of the Universe. But these questions should be considered in further studies. 

We have also addressed the issue of how the effects related to the nonminimal coupling could be manifested in 
laboratory experiments. One such effect is the occurrence of extra accelerations whose values must lie within the 
restrictions imposed by various experiments. Proceeding from a consideration of the equation of motion for a fluid 
element in a given gravitational field, in the Appendix the restriction imposed by the Eotvos-type experiments is 
derived [see Eq. (IA.5I) ]. Despite the fact that the nonminimal function / obtained here leads to the unusual behavior 
of the characteristics of compact configurations, we use this function to demonstrate the theoretical possibility of 
how one can estimate extra accelerations which appear in the Earth's atmosphere because of the presence of the 
nonminimal coupling. It is shown that the condition (|A.5[) is satisfied in principle [see Eq. (IA.14|) ] . In future work, 
we plan to perform a consideration of these questions using such a nonminimal coupling function, which gives a more 
adequate description of the characteristics of compact configurations. 



The author is grateful to the Research Group Linkage Programme of the Alexander von Humboldt Foundation for 
the support of this research. This work was partially supported by a grant in fundamental research in natural sciences 
by the Ministry of Education and Science of Kazakhstan. 



In this appendix we address the question of how the effects related to the presence of the nonminimal coupling in 
the system described by the Lagrangian ([T]) could manifest themselves in laboratory experiments. 

To illustrate this, consider the equation of motion for a fluid element in a gravitational field. We derive this equation 
using the energy-momentum tensor ([3]) which can be rewritten in the form 
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Appendix A: Effects of the nonminimal coupling and laboratory constraints 




(A.l) 



where 



= (e+ P ) 
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Here the total energy-momentum tensor is divided into two parts denoted by the indexes (f) and (sf) which refer to 
the fluid and the scalar field, respectively. On operating a covariant derivative on (|A.1I) , we have 

By virtue of the scalar field equation ([5]), the third term on the right-hand side of this expression vanishes. Projecting 
two remaining terms onto the direction of the fluid 4- velocity it M , one can obtain the following equation of motion for 
a fluid element: 

dV< dx^dx^ _ p ; „+p(ln/) ; „ 

ds* +L °u ds ds - £+p ^ « 9 )■ (A.2j 

The first term in the numerator on the right-hand side, p. v , is the pressure gradient that also appears in usual 
non-chameleon models and describes the force exerted on a fluid element due to the fluid pressure (it is not to be 
attributed to the coupling between the matter and the scalar field, nor does it signal any new effect). The presence 
of this term forces particles of matter, within which there are pressures that generally move along lines that are not 
geodesies. The second term in the numerator, p(hxf) ;u , is provided by the presence of the nonminimal coupling in 
the system. One can see that its influence on the motion of a test particle is determined both by the magnitude of 
the fluid pressure and the gradient of the nonminimal coupling function /. In the general case when magnitudes of 
the pressure and the gradient of / are not small, one might expect that substantial extra accelerations in the system 
described by the Lagrangian (|TJ) will appear. 

The presence of such extra accelerations implies that the theory may exhibit violations of the weak equivalence 
principle. The most stringent constraints on possible violations of the equivalence principle derive from experiments 
performed on Earth. These experiments test the universality of free fall, the prediction that all bodies in a uniform 
gravitational field have exactly the same gravitational acceleration. The universality of free fall has been tested to 
roughly 1 part in 10 12 26] using laboratory test bodies in the gravitational fields of the Earth and the Sun. Another 
sort of experiments is connected with precise measurements by using lunar laser-ranging to compare the accelerations 
of the Earth and Moon toward the Sun. These experiments constraint the difference in free-fall acceleration to be less 
than approximately one part in 10 13 (see, e.g., Ref. [11]). 

Then, proceeding from the requirement to satisfy the aforementioned laboratory experiments, one must make sure 
that the effects related to the presence of the nonminimal coupling do not lead to a contradiction with the experiments. 
To check this, we will estimate below the value of extra acceleration which appears in the presence of the nonminimal 
coupling with the function / given by the expression (|2ip. But before doing that, it will be useful to notice that 
since we consider here the case of low-pressure approximation, when f(x) ~ f c — const, the gradient of the function 
/ is virtually equal to zero along the radius of the configuration. This allows us to expect that the influence of the 
term p(ln/) ; „ in Eq. (|A.2j) will be small compared with the gradient term Indeed, the numerical calculations 
show that the smaller the value of the relativistic parameter er s , the smaller the value of the term p (ln/) ; „ compared 
with p. v . One would expect that in the nonrclativistic limit, when <j s — > 0, the extra acceleration associated with the 
nonminimal coupling will be negligibly small compared with the acceleration provided by a gravitational field. 

It is natural to assume that the nonrelativistic approximation used in the present paper is also relevant in describing 
relatively small objects such as planets. Consistent with this, let us estimate the value of extra acceleration in a 
nonrelativistic ("Newtonian") system in the presence of a nonminimal coupling such as pf(ip). In doing so, we will 
compare the gravitational acceleration of a fluid element with the extra acceleration associated with the nonminimal 
coupling. To do this, let us find from Eq. (IA.2I) the expression for the usual three-acceleration in the nonrelativistic 
approximation p<£ 

atot « a N ( 1 + — ] , (A. 3) 

\ a N J 

where a to t is the total radial acceleration defined as a sum of the Newtonian acceleration, — —dip/dr (ip is the 
Newtonian gravitational potential), and the extra acceleration provided by the chameleon, a c h = —(jp/p){d\xif/dr). 
To compare the expression (|A.3[) with experimental data, it is convenient to use the Eotvos parameter, 

_ Joi -02I 

V = 2, ; r, 

|oi + a 2 | 

which describes the difference in relative free-fall acceleration for two bodies of different composition having the 
accelerations ai and 0,2- As mentioned above, r\ is experimentally restricted to be less than 10 -13 . Taking into 
account the smallness of 77, one can find from the above expression: 

a x « 02(1 + 77). (A.4) 
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Compare the expression (IA.3I) with (IA.4I) : whereas the former gives the magnitude of the total acceleration of a fluid 
element in the presence of the nonminimal coupling, the latter describes the difference between the accelerations of two 
test bodies which are manufactured to have different compositions. But in both cases the difference in accelerations 
is estimated, related to either the nonminimal coupling (in the first case) or a different composition of bodies (in the 
second case). It is clear that independently of the fact which case is under consideration, the experiment requires 
that the difference in accelerations be less than r\ rj 1CP 13 . Then, by identifying formally atot = ffli, ajv = ^2, and 
ich/ajv = we have the following restriction on the extra acceleration from (|A.3[) : 

p d In f „ 1 a , . , 

a ch = — -L < 10- 13 a N . A.5 
p dr 

Here the pressure p and the mass density p refer to a fluid element falling onto an attractor. The interaction of this 
matter with a chameleon field through the function / results in the extra acceleration a c h, the value of which must 
satisfy the condition (|A.5[) to be consistent with the experiment. 

It is seen that the condition (|A.5[) is trivially satisfied if a continuous medium, located in the Earth's gravitational 
field, is pressureless. Such a case corresponds, for instance, to dustlike matter. If the pressure is nonzero, the problem 
requires a separate analysis for each particular case. Then, from the point of view of experiment, each element of 
such a fluid can be considered as a test particle moving in a given gravitational field. In this case the magnitude of 
extra acceleration will depend on the composition of this medium, which, in turn, is determined by use of various 
equations of state p — p(p) for different forms of matter. Then, aside from the fact that the extra acceleration does 
take place for one type of matter, there will also be additional differences in accelerations for fluids having different 
compositions. 

1. Extra accelerations in the Earth's atmosphere 

To illustrate the above considerations, let us now demonstrate how one can estimate the value of extra accelerations, 
which could appear, for instance, in the Earth's atmosphere. As is known, the Earth's atmosphere is not homogeneous 
and the air density decreases with altitude; the dependence of the density p on altitude R can be approximated by 
the barometric formula 

P = p e-( R - R ^l A . (A.6) 

Here po is the density at sea level, R = Rq, and A is the so-called scale height of the standard atmosphere, which 
at the Earth's surface is approximately equal to 8.5 km. For our purposes, we need only to pick this layer of the 
atmosphere, i.e., in further calculations we assume that i?0 < R < (-R© + A). This choice of p implies that the air 
temperature does not change with altitude. This allows one to describe the atmosphere by using the equation of state 
for an isothermal fluid given in the form 

p = Kp (A.7) 

with K = po/ po = c„, where c a is the speed of sound in the air, and po is the pressure at sea level. Using this equation 
of state and introducing the dimensionless variables 



e"* = A C=£ *(0 = ^P*(r), (A.8) 
pa L K 

where L = \J K / {AnG po) has dimensions of length, one can obtain from the Lagrangian ([1} the following nonrelativistic 
equations describing the isothermal fluid in the presence of a chameleon scalar field [19| : 

1 d ( ' 2 dx 



c^cv c ^r /e " x ' (A - 9) 

The first equation is the equation of hydrostatic equilibrium of a fluid in a Newtonian gravitational field. Substituting 
the given distribution of density (|A.6[) in (|A.9[) . one can find 

r , „ exp [8(( - C ffi )l 
/[0(C)] =2/3 ^® n , 
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where j3 — L/A ss 10 3 . Using this expression in (|A.10[) and taking into account that j3 > 1, one can obtain an 
equation for the scalar field in the form 



2fP_ 



2 vt , C r 

where the prime denotes differentiation with respect to £. The first integral of this equation is 



6* = ^ 
C 



The value of the integration constant C\ is determined from the condition that at the edge of the atmosphere 4>' ~ 
(since outside the atmosphere the scalar field makes a smooth transition into the cosmological chameleon field). 
Representing £ in a form £ = £0 +6, where 5 = A/L « 10~ 3 , and substituting this in the above expression, we finally 
get 



c 



4(5 \ 

1 + 



(A.ll) 



The solution of this equation gives the distribution of the scalar field in the Earth's atmosphere. We will seek a 
solution starting from the Earth's surface, £ = (^g, and choosing the initial value of 4>(C(&) i n such a way as to get 
approximately the cosmological value of the field at the outer boundary of the atmosphere. 

Since the cosmological scalar field, <^>(£), and that used here, </>(C): are measured in different units, then, to find the 
relation between them, compare their definitions given by the expressions (fTT]) and (|A.8|) , respectively. Thus we have 



K 



0(O = V2^0(C). 



(A.12) 



Substituting this expression in (|2"Tj) . we get 



Oa [ er(ri + 1) 
— O 



-2(n+l) 



(A.13) 



Next, the experimental restriction from (|A.5[) on the extra acceleration near the Earth's surface can be rewritten, 
by using the dimensionless variables from (|A.8|) . as follows: 



din/ 



< 10 



-10 



(A.14) 



Here we have used the following numerical values: a at « 10 m s -2 and K — po/po ~ 10 5 m 2 s~ 2 . Then, evaluating </)(() 
numerically by using Eq. (|A.11|) . and substituting this in (|A.13p . it must be checked whether the condition (|A.14j) is 
satisfied or not. The calculations indicate that the value of the logarithmic derivative of (|A.13|) . 



din/ 



1)^^(0 coth 



near the Earth's surface is of order 10~ 10 for the values of n used in the present paper. Thus, the condition (|A.14[) 
may be considered as satisfied, in principle. The estimate for the condition (|A.5[) for other forms of matter requires 
a separate analysis. 
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